Spectral analysis of gluonic pole matrix elements for fragmentation 
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The non- vanishing of gluonic pole matrix elements can explain the appearance of single spin asym- 
metries in high-energy scattering processes. We use a spectator framework approach to investigate 
the spectral properties of quark-quark-gluon correlators and use this to study gluonic pole matrix 
elements. Such matrix elements appear in principle both for distribution functions such as the Sivers 
function and fragmentation functions such as the Collins function. We find that for a large class 
of spectator models, the contribution of the gluonic pole matrix element in fragmentation functions 
vanishes. This outcome is important in the study of universality for fragmentation functions and 
confirms findings using a different approach. 



I. INTRODUCTION 

In high-energy scattering processes the structure of 
hadrons is accounted for using quark and gluon corre- 
lators, forward matrix elements of non-local quark and 
gluon operators between hadronic states. Making an ex- 
pansion in the (inverse) hard scale, the relevant compo- 
nent of the momentum of partons (quarks and gluons) 
is the one collinear to hadrons and correspondingly the 
non-locality in the matrix elements is restricted to the 
light-cone. Moreover, all leading dynamical effects come 
from two-field configurations at two light-like separated 
points, which are easily interpreted as parton densities 
or parton decay functions [Ullj]- These are the parton 
distribution functions depending on the momentum frac- 
tion x relating the parton momentum k — x P to the 
hadron momentum P or the fragmentation functions of 
partons into hadrons depending on the momentum frac- 
tion z, relating the parton momentum k and the hadron 
momentum P = z k. 

At sub-leading order in the hard scale or when ex- 
plicitly measuring transverse momenta, other matrix el- 
ements become important such as the three-parton cor- 
relators containing parton fields at three different space- 
time points with light-like separations and two-parton 
correlators with also transverse separation (light-front 
correlations). These latter (light-front) correlators are 
described in terms of transverse momentum dependent 
(TMD) distribution and fragmentation functions, which 
are sensitive to the intrinsic transverse momenta of par- 
tons in hadrons, k = x P + k T in a frame in which the 
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FIG. 1: The graphical representation of the quark-quark cor- 
relator in the case of distributions of partons with momentum 
k in a hadron with momentum P (a) and the spectator model 
description (b). 



hadron does not have transverse momentum (P T = 0) or 
for fragmentation k = - P + k T . In this case one often 
refers to the hadron transverse momentum P± = —zk T 
(in a frame in which the parton does not have a trans- 
verse momentum (kj_ = 0)). 

In this paper, we will investigate multi-parton cor- 
relators with one additional gluon in which the zero- 
momentum limit will be studied 0, i, i, H, 0, H, H, E31 • 
These are so-called gluonic pole matrix elements or 
Qiu-Sterman matrix elements, that have opposite time- 
reversal (T) behavior as compared to the matrix elements 
without the gluon. Such matrix elements involving time- 
reversal odd (T-odd) operator combinations are of inter- 
est because they are essential for understanding single 
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FIG. 2: The graphical representation of the correlator in the 
case of fragmentation of partons with momentum k into a 
hadron with momentum P (a) and the spectator model de- 
scription (b). 



spin asymmetries at high energies. In order to under- 
stand the basic features of these matrix elements we per- 
form a spectral analysis by modeling the distribution and 
fragmentation functions under reasonable assumptions. 

For the correlators, depicted in Figs. QJa) an d H^a), 
one has expressions in terms of matrix elements of bilo- 
cal operators that are frequently used as a starting point 
in modeling distribution and fragmentation functions. 
In particular the spectator model, pictorially given in 
Figs, [lib) and[^b), has become fairly popular, because 
it is easy, flexible and intuitively attractive. On the other 
hand, one should be very careful, because the predic- 
tive power depends on limiting oneself in the choice of 
spectator (e.g. a diquark with fixed mass in case of the 
nucleon) and using simple vertices. In fact making a 
spectral analysis of the spectator and allowing for the 
most general vertices one would lose all predictive power. 
Having said these words of caution, we will investigate in 
this paper differences between distribution and fragmen- 
tation functions using a spectral analysis and using phys- 
ical intuition in restricting the momentum dependence 
and asymptotic behavior of the vertices. In this context, 
the relevant gluonic pole matrix elements that we want 
to study are &c{k,k — k\) and Aa(k,k — k\) shown in 
Figs. [3] and 2] Of these matrix elements only the de- 
pendence on the collincar components x and x\ in the 
expansion of the momenta is needed. We find that while 
both §q(x,x — x\) and Aq(x,x — x\) are nonzero, tak- 
ing the limit x\ — > x, < I ) g(x,x) remains non-zero, while 
Aq(x,x) vanishes. 

The vanishing of the T-odd gluonic pole matrix ele- 
ments is important in the study of universality of TMD 
distribution functions (DFs) and fragmentation functions 
(FFs). For the collinear case T-symmetry can be used as 



a constraint on the parton correlators, limiting the DFs 
to T-even ones. Such a constraint does not apply for the 
fragmentation correlator because the final state hadron 
is part of a jet and as such not a plane wave, allowing 
both T-even and T-odd FFs. But for spin and spin 
1/2 hadrons no T-odd functions appear at leading twist 
(leading order in the hard scale) [Tllll^]. Including trans- 
verse momentum dependence, both the distribution and 
fragmentation correlators ($ and A) are no longer con- 
strained by T-symmetry. The reason is that the appro- 
priate color gauge invariant operators in the correlator, 
in particular the gauge links, are not T- invariant. For 
DFs this provides a mechanism leading to T-odd func- 
tions, such as the Sivers function [3]. But, for FFs there 
are now in principle two mechanisms leading to T-odd 
functions [l4j. 

A nice feature, however, is that the two mechanisms 
leading to T-odd functions can be distinguished. The T- 
odd operator structure can be traced back to the color 
gauge link that necessarily appears in correlators to ren- 
der them color gauge-invariant. But the operator struc- 
ture of the correlator is also a consequence of the nec- 
essary resummation of all contributions that arise from 
collinear gluon polarizations, i.e. those along the hadron 
momentum. How this resummation takes effect is a mat- 
ter of calculation. The result is a process dependence in 
the path in the gauge link. After azimuthal weighting of 
cross sections one simply finds that the T-odd features 
originating from the gauge link lead to specific factors 
with which the T-odd functions appear in observables. 
Comparing T-odd effects in DFs in semi-inclusive deep 
inelastic scattering (SIDIS) and the Drell-Yan process 
one finds a relative minus sign [TEI ] . Similarly, comparing 
T-odd effects in FFs in SIDIS and electron-positron an- 
nihilation one also finds a relative minus sign, at least for 
the T-odd effect originating from the operator structure 
(gauge link) The effect coming from the hadron-jet 
final state not being a plane wave will not lead to process 
dependent factors. 

Before coming to the significance of our work, we need 
one more ingredient. We already mentioned that relating 
T-odd effects in different processes, requires azimuthal 
weighting, which projects out the transverse momentum 
weighted parts of the correlators $ and A, referred to as 
transverse moments $g and Ag, respectively. The T-odd 
operator parts are precisely the soft limits (k\ — > 0) of 
the gluonic pole matrix elements [l^|. Thus 

$a = $a + 7r$ G (fci = 0), (1) 
A s = A a + ^A G (fc 1 = 0), (2) 

When Ac(ki — 0) is zero, there still are T-odd FFs con- 
tained in Ag. They appear in the matrix elements of the 
T-even operator combination in Ag involving a hadron- 
jet state (non-plane- wave) and they are process indepen- 
dent, for instance the T-odd Collins function In con- 
trast T-odd DFs in $a only can come from $G(ki = 0). 
These DFs can still be universal but appear with calcu- 
lable process dependent gluonic pole factors [13, [3 ■ 
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FIG. 3: The graphical representation of the quark-quark- 
gluon correlator <&g for the case of distributions including a 
gluon with momentum ki (a), and the possible intermediate 
states (b) and (c) in a spectator model description. 



The above had already been shown in model calcu- 
lations for SIDIS and e + e _ annihilation I19l [2011 and 
more recently in hadron-hadron scattering [2X| . I'2'2| |. In 
these calculations the authors look at the full process 
in the model, carefully studying the cuts, rather than 
concentrating on the soft part only. By contrast, we 
look at the soft part only, that is, the multi-parton cor- 
relators $g and Aq- In principle such an approach 
of only looking at the soft part is also possible start- 
ing with TMD dependent two-parton correlators <E> and 
A [H, [H [H H [H [H, HI, M M Hi- To generate 
T-odd contributions for TMD functions in such an ap- 
proach one needs to go at least to one-loop [H[ perform- 
ing a perturbative expansion on the gauge link 0, [25| . 
Extraction of the complete T-odd gluonic pole contribu- 
tion is difficult. Our approach, starting directly with 
the color gauge invariant multi-parton correlator (that is 
having re-summed the gauge link [TJ]) has the advantage 
we can work with tree-level matrix elements and just per- 
form a spectral analysis to extract the T-odd gluonic pole 
contributions. 

In the next section we give details on the gluonic pole 
matrix elements, followed by the spectator model ap- 
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FIG. 4: The graphical representation of the quark-quark- 
gluon correlator Ag in the case of fragmentation including 
a gluon with momentum ki (a) and the possible intermediate 
states (b) and (c) in a spectator model description. 



proach. In the concluding section, we outline the short- 
comings of our approach and also discuss possibilities to 
use the approach for more detailed estimates for the T- 
odd FFs. 

II. GLUONIC POLE MATRIX ELEMENTS 

At high energies, it is useful to make a Sudakov de- 
composition of the momenta of active partons, k = 
x P + a n + k T . The Sudakov vector n is an arbitrary 
light-like four-vector n 2 = that has non-zero overlap 
P ■ n with the hadron's momentum P. We will simply 
choose P ■ n = 1. In a hard process, the Sudakov vec- 
tor incorporates the presence of other momenta that are 
hard with respect to the hadron under consideration, e.g. 
n w P'/P ■ P'. We can now also work with light-cone co- 
ordinates. Including mass effects one would have n_ = n 
and n+ = P 



i M 2 n and with k* 1 



k ■ 7i zp we have 



k+ 
k~ 



k ■ n - 
k-P 



\xM 2 



\xM 2 . 



(3) 
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Vectors in the transverse plane can be obtained from the 
transverse projector, g? v — g^ v — n^n"} . 

As the effects of the component fc _ will appear sup- 
pressed by two powers of the hard scale as compared to 
the collinear term, it is integrated over and one considers 
quark-quark correlators on the light-front (LF: £ • n = 0) 



(2tt) 3 

x(P\^(0)U m MO\P)\ LF ■ (5) 

The Wilson line or gauge link Wu.g = 
■Pexp[— igf c ds-A a (s) t a ] is a path-ordered expo- 
nential along the integration path C with endpoints at 
t] and £. Its presence in the hadronic matrix element is 
required by gauge-invariance. In the TMD correlator ([5]) 
the integration path C in the gauge link is process- 
dependent. In the diagrammatic approach the Wilson 
lines arise by resumming all gluon interactions between 
the soft and the hard partonic parts of the hadronic 
process [H, H, M, HI • 

Collinear quark distribution functions are obtained 
from the TMD correlator after integration over p T , 

$(s) - J d 2 k T $M(x,k T ) 
'd(ti-P) 



2tt 



e*^- p (P\mU(o^(t)\P)\ LC ■ 

(6) 



The non-locality is restricted to the light-cone (LC: £-n — 
£ T = 0) and the gauge link is unique, being the straight- 
line path along n. In azimuthal asymmetries one needs 



the transverse moments contained in the correlator 

$a M 0) = J d2k T K ^ [U] {x,k T ) . (7) 

The TMD correlator, expanded in distribution functions 
depending on x and k\ contains T-even and T-odd func- 
tions, since the correlator is not T- invariant, which is 
attributed to the gauge link that depending on the pro- 
cess, accounts for specific initial and/or final state inter- 
actions depending on the color flow in the process. For 
the collinear case, the link structure becomes unique in 
the case of integration over k T (Eq. [5]). For spin and 
spin 1/2 the quark and gluon correlators that appear at 
leading order in high energy processes contain only T- 
even operator combinations. Evaluated between plane 
waves one only finds T-even functions depending on x in 
the parameterization. 

For the collinear weighted case, the transverse mo- 
ments in Eq. Q one retains a nontrivial link-dependence 
that prohibits the use of T-invariance as a constraint. It 
is possible, however, to decompose the weighted quark 
(and also gluon) correlators as 



(8) 



with calculable process-dependent gluonic pole factors 
C [ q ] and process (link] independent correlators $a and 
<&G. The correlator &g contains the T-even operator 
combination, while $g contains the T-odd operator com- 
bination. The latter is precisely the soft limit x\ — > of 
a quark-gluon correlator &a{x, %i) of the type 



2vr 2tt 

The universal T-odd distribution functions in the parameterization of &g(%i %) appear in T-odd observables such as 
single spin asymmetries with the specific gluonic pole factors from Eq. [8] 

The situation for fragmentation functions is different. The TMD fragmentation correlator depending on the collinear 
and transverse components of the quark momentum, k — - P + k T + crn, is given by 14 1 

A^(z,M = Y,J d ^' { 2^ T e ik <{0\U m MO\P,X){P,X\^(0)\0)\ LF . (10) 

The collinear, A: T -integrated correlator 



A(z) = f d 2 k T AM(z,k T )=J2 I e lz ' 1 ^ p ^0\U^ ] MO\PX)(pX\^(0)\0)\ LC , 
J x J K 



ill) 



only contains a T-even operator combination. Nevertheless one could in principle have T-even and T-odd fragmen- 
tation functions depending on z since the hadronic state \P,X) is an out-state, which is not T-invariant. For spin 
and spin 1/2 hadrons no T-odd function appear at leading twist because of other constraints. At sub-leading twist 
they do appear (37j |. 

In the transverse moments obtained after fc T -weighting, 



A a 9 M(z) = J d 2 fc T fc?A^(z,fc r ) = Ag(i) +cM*Aa(i,§) 



(12) 
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the two link independent correlators Ag and Aq contain again a T-even and T-odd operator combination, respectively. 
The gluonic pole correlator is again the soft limit, z^ 1 = x% — ► 0, of the quark-gluon correlator 



A' 



LC 

(13) 



Because of the appearance of hadronic states \P, X), each 
of correlators in Eq. [T2] contains in principle T-even and 
T-odd functions. Rather than having a doubling of T- 
odd functions, we will show in a spectator model ap- 
proach that Aq(x,x) = 0, which implies that T-odd 
fragmentation functions in the transverse moments only 
come from Ag, which appear with a universal strength 
(no gluonic pole factors). We will show this in a spec- 
tator approach starting with the collinear quark-gluon 
correlators in E qs. El and I13|r ather than the model ap- 
proaches m, m, (Mil, ejee m m m hi that looked 

at the transverse momentum dependent quark correlators 
in Eqs and [TO] 



III. THE SPECTATOR MODEL APPROACH 

In a typical spectator model approach to distribution 
or fragmentation correlators one considers a spectator 
with mass M s . The result for the cut, but untruncated, 
diagrams, such as in Figs. [T] and [2] are of the form 



$(x,fc T )~ / d(k-P) 



F(k 2 ,k- P) 



(fc 2 



+ ie) 



rS ((fc - Pf - Ml) , 
(14) 



where F(k 2 , k-P) contains the numerators of propagators 
and/or traces of them in the presence of Dirac Gamma 
matrices, as well as the vertex form factors (see for ex- 
ample [38]). The explicit momenta, using light-cone co- 
ordinates [p~ , p + , p T ] as discussed in the beginning of the 
previous section, are 



P = 
P-k 

k = 



M 2 

— 1,0, 
M 2 - 



1 



2(1 -x) '" 
\(l-x)M 2 - M 2 



2(1 -x) 



(15) 
(16) 

(17) 



In the above the delta function constraint in Eq. [TJ] 
has been implemented. One finds that the numerator 
F(fc 2 , k-P) = F(x, fc 2 ,) and hence 



$(x, kj 



(1 - a;) 2 F(a;,fc r ) 
(H 2 (x)-k 2 f ^ 



(18) 



ith 



H 2 {x) = xM 2 s + (1 - x)m 2 - x(l - x)M 2 



(19) 



Note that fc 2 = — fc^ < 0. The details of the numerator 
function depend on the details of the model, including the 
vertices, polarization sums, etc. These must be chosen in 
such a way as to not produce unphysical effects, such as 
a decaying proton if M > m + M s , thus m in Eq. [Ml must 
represent some constituent mass in the quark propagator, 
rather than the bare mass. The useful feature of the 
result in Eq.[T5]is its ability to produce reasonable valence 
and even sea quark distributions using the freedom in 
the model connected to an intuitive picture. The results 
for the fragmentation function in the spectator model is 
identical upon the substitution of x = l/z. 

Next we turn to the same spectral analysis of the glu- 
onic pole correlator using the picture given in Figs. [3] for 
distribution functions and the picture given in Figs. [H 
for fragmentation functions. Again, we only need to in- 
vestigate one of the cases. We parameterize the gluon 
momentum as 



fci = 



fc, , x\ , k\i 



(20) 



where ky = k% -P— h x\ M 2 will be the first component to 
be integrated over. The relevant momenta (implementing 
the on-shell condition for P — fc) are 



fc — fci = 
P-k + kx = 

P-k! = 



(1 - x) M 2 - M 2 + fc 2 
fci + „/-, s i x x\,k T ki 2 

2(1 — x) 

• _ M 2 -fc 2 1 

fc l + TTT, TV i i- — X + X\, —k T + fcii 



2(1 -x) 



M 2 



1 - x%, -ki T 



(21) 
(22) 
(23) 
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The basic result for the quark-gluon correlators &g( x i % ~ %i, k T , k T — k\ T ) becomes 



dk. 



F 1 (fc-L ,x,xi,k T ,k lT ) 



(k 2 -m 2 )\J 2ni (fc 2 - ml + ie)((k - ki) 2 - m 2 + ie)((P - k + ki) 2 - + ie) 



dk 1 



F 2 (k x ,X,Xx,k T ,k lT ) 



1 - x 



2iri (fef - mf + ie)((k - h) 2 - m 2 + ie)((P - h) 2 - M 2 2 + ie) j 

dki F l (fcf , x, Xx, fc T , fc lT J 

{n 2 -k 2 )\J 2vr i (xi - A\ + ie)((xi - x) fcf - A 2 + ie)((l - x + xi) fcf - S x + ie) 

, X, Xi, fcy, fc 2 T ) 

27ri (xi —Ai + ie)((xi - x) fcf - A 2 + ie)((xi - 1) fcf - S 2 + ie) 

I 



(24) 



where, as before, F{ {k 1 , y, Xi , /c ^., /civ) contain numera- 
tors and vertex functions [24, 25}. We use the quantities, 



2 A\ = ml - k 



It 1 



(25) 



2A 2 = m A - (x-xx)M* + 



2 , X X l 



1 -X 



(M 2 



\k T -k lT ) 2 , (26) 



2 Si = M. 



.si 



(M s 2 - fc 2 



(fe T - fc lr ) 2 , (27) 



1 -x 

2 ,i „ \ ^2 ,.:2 



2S 2 = Mi 2 -{l- xx)M z -kf T 



(28) 



These quantities depend on spectator masses, momentum 
fractions (x and Xi) and the transverse momenta (k T 
and fci r ). Besides the spectator mass M s , two additional 
spectator masses M s \ and M s2 appear (see Figs.lSJb) and 
(c) or Figs. 0Jb) and (c) and compare with the starting 
expression in Eq. I24[) . The quantity /i 2 is the same one as 
given in Eq. 1191 Assuming that the numerator does not 
grow with k^ one can easily perform the fc-j" integrations. 
We will for simplicity assume that the Fi are independent 
of k x , in which case we obtain for x > x\ > (see 
Appendix) , 



&g{x, x — X\) = I d 2 k T d 2 ki T 



(1 -x)i 7, i(x,xi,fc r ,fci T ) 



(/i 2 - k 2 ) (x A 1 + Xl (A 2 - A,)) (x x (S x - A 2 ) - xSi - (1 - x) A 2 ) ((1 -x)A 1 + x x (A, - B{)) 

(xi (Si - A 2 ) - x Si - (1 - ar) A 2 ) xi<9(xi) 

+ ((1 -x)Ai + x x (Ax - Si)) (xi - x)6(xi - x) 

+ (xA x + Xl (A 2 - A x )) (1 - x + x x )0(l - x + Xl ) 

(1 - x)F 2 (x,xx,k T ,k lT ) 

( M 2 - k 2 ) (x A x + Xl {A 2 - Ax)) (A 2 -xB 2 + x x (S 2 - A 2 )) (xi {Ax - B 2 ) - Ax) 

(A 2 -xB 2 + Xl (S 2 - A 2 )) xxO(xx) 

+ (xi (Ax - S 2 ) - Ax) (x x ~ x)6(xx - x) 

+ (xAx+xx{A 2 -Ax)){xx-l)6{xx-l) j (29) 

Taking the limit xi — * we get the gluonic pole correlators, for distribution functions (0 < x < 1), 

(1 - x)Fx(x,O,k T ,kx T )0(l -x) 



$g(x, x) 



d 2 k T d 2 kx T 



(n 2 -k 2 )(xBx + (l-x) A 2 ) Ax 



(30) 
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and for fragmentation functions (x = X/z > 1) 

A G (x,x) = 0. (31) 

We note that this result depends on the assumption that 
the numerator does not grow with . If this is the case 
one would find integrals of the type Q\n (also given in 
the Appendix) rather than those of the type 0° n and one 
does not get the required x\ 9(xi) behavior in the calcu- 
lation. In models, terms proportional to fcf ~ fci • P 
may easily arise from numerators of fermionic propaga- 
tors [H, [39[ which in turn may easily be suppressed by 
form factors at the vertices. To prove a proper behavior 
within QCD one would need to study the fully unintc- 
grated correlators such as e.g. in Ref. [4(| and show that 
they fall off sufficiently fast as a function of k\ ■ P. 

IV. CONCLUSIONS 

In this work, we have investigated the gluonic pole 
contributions to the distribution and fragmentation func- 
tions. Instead of doing a quantitative analysis involving 
details of a phenomenological model, we limit ourselves 
to a spectral analysis within the spectator framework, in 
order to understand the basic features of these quantities. 
The advantage is that we are able to investigate only the 
soft parts at tree level and take the zero momentum limit 
of the gluon involved. We simply assume that masses and 
vertices do not spoil our analysis, which implies limits 
on the mass distributions of the spectators, use of ver- 
tices that cancel the bare-mass poles in the quark and 
gluon propagators and behavior of vertices that assures 
sufficient convergence of integrations. We find that un- 
der realistic assumptions, the gluonic pole contributions 

I 



for fragmentation correlators vanish whereas these con- 
tributions do not vanish for distribution correlators. The 
result for fragmentation correlators at nonzero gluon mo- 
mentum is nonzero. We stress that this certainly is not 
yet the full proof that gluonic pole matrix elements van- 
ish in the case of fragmentation. However, we consider 
this analysis as a step towards such a proof and the pos- 
sible direction to obtain such a proof by only consider- 
ing the appropriate color gauge-invariant soft matrix ele- 
ments. Such a proof is important as it eliminates a whole 
class of matrix elements parameterized in terms of T-odd 
fragmentation functions besides the T-odd fragmentation 
functions in the parameterization of the two-parton cor- 
relators. For instance, only one of the contributions to 
the spin asymmetries considered for jet-hyperon produc- 
tion in Ref. [4l| remains. Moreover, the remaining frag- 
mentation functions appear with the standard partonic 
cross section, so no gluonic pole cross sections need to be 
considered for fragmenting final-state partons, limiting 
these considerations to the distribution functions involv- 
ing initial-state partons. 

As mentioned, the results in this paper may point the 
way to find a full proof of Aq(x, x) — 0. The approach 
taken here to look at tree- level three-parton correlators, 
also can be used for explicit model calculations for T- 
odd distribution functions originating from gluonic pole 
matrix elements and the investigation of their effects in 
single spin asymmetries. 
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APPENDIX A: USEFUL INTEGRALS 



Often it is useful to attack integrals containing propagator poles via light-cone variables, leading to integrals of the 
type 



+ ie) ni (ax 2 - 1 + ie)" 2 
for which easy reduction rules exist 42] . We need specifically 



(Al) 



Q u {xi,x 2 ) = = , (A2) 

X\ — x 2 xi- x 2 

Q°m(xi,x 2 ,x 3 ) = -. X2 r e° n (x 2 ,x 3 ) - Xl , eUxi,x 3 ) (A3) 

(Xl-X 2 ) (X1-X2) 

(x 2 - x 3 ) x x 6{x x ) + (x 3 - xi) x 2 0(x 2 ) + (xi - x 2 )x 3 8(x 3 ) 



(xi - x 2 )(x 2 - x 3 )(x 3 - Xl) 
(x 2 - x 3 ) 6{xi) + (x 3 - xi) 6(x 2 ) + (xi - x 2 ) 9(x 3 ) 



©inOi, x 2 , x 3 ) = - [Q n(x 2 ,x 3 ) - 9°! (£1,2:3)] (A5) 



(xi ~ x 2 )(x 2 - x 3 )(x 3 - Xl) 



(A6) 



Including arbitrary pole positions, one finds 
da 1 



2iri (axi — A\ + ie)(ax 2 — A 2 + ie) 



x\ - x 2 



e° n (x u x 2 ) 



x\A 2 - x 2 Ai 
0(xi) - 6(x 2 ) _ 6{ Xl ) 9(-x 2 ) - 0(- Xl ) 0(x 2 ) 



xiA 2 - x 2 Ai 



xiA 2 - x 2 Ai 



1 



da 

2ir i (ax\ — A\ + ie) (ax 2 — A 2 + ie) (axi — A3 + ie) 



1 



XiA 2 - x 2 Ax 



x 2 



x 2 - x 3 
x 2 A 3 - x 3 A 2 



Q° 11 (x 2 ,x 3 ) - x\ 



xi - %3 
x ± A 3 - x 3 Ax 



9?! (0:1,0:3) 



_ (x 2 A 3 - x 3 A 2 ) xi 0(xi) + {x 3 Ai - X1A3) x 2 9(x 2 ) + (x 1 A 2 - x 2 A x ) x 3 6(x 3 ) 
{x x A 2 - x 2 A 1 )(x 2 A 3 - x 3 A 2 ){x 3 A x - x x A 3 ) 
da a 

2ir i (ax x — A x + ie) (ax 2 — A 2 + ie) {ax x — A 3 + ie) 

x 2 - x 3 x x - x 3 
A 2 3 — 611(0:2,0:3) - A x — t 811(0:1, 0:3) 



2:1^2 - 0:2^1 



x 2 A 3 - x 3 A 2 



xiA 3 - x 3 Ai 



(x 2 A 3 - x 3 A 2 ) 6{x t ) + {x 3 A x - Xl A 3 ) 9(x 2 ) + ( Xl A 2 - x^) 9{x 3 ) 
{x x A 2 - x 2 Ai){x 2 A 3 - x 3 A 2 ){x 3 Ai - x x A 3 ) 

I 



(A7) 
(A8) 

(A9) 
(A10) 

(AH) 
(A12) 
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